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Introduction. Properties of isentropic ex- 
pansion of warm dense matter (WDM), which 
is created by intense energy flux (strong shock 
compression pQ or intence isochoric heating by 
laser or/and heavy ion beam [2]) are often 
discussed in situation when: (i) - thermody- 
namic trajectory of such expansion crosses bin- 
odal of liquid-gas phase transition (i.e. boiling 
curve), and (it) - isentropic expansion in two- 
phase region is going along equilibrium branch 
(not metastable one) of two-phase mixture isen- 
trope. Possibility of formation of remar cable 
"plateau" in density and pressure profiles was 
studied previously in [3 (and claimed in [4]) 
where the problem of metal surface ablation un- 
der the action of strong laser radiation was con- 
sidered. Similar plateau was announced also in 
hydrodynamic symulations of equilibrium isen- 
tropic expansion with real two-phase EOS of 
some metals [5]. 

Peculiarity of "freezing" of finite portion of 
expanding matterial in extended and uniform 
state ("boiling" layer) corresponding to binodal 
of gas-liquid or /and other phase transitions, - 
so called "phase freezeout" - and prospects of 
applications of this phenomenon for intended 
generation of uniform and extended zone of pre- 
viously unexplored states of matter just at the 
two-phase bundry were discussed in [6]. 

Features and parameters of such "boiling" 
liquid layer, which could arise under isentropic 
expansion of WDM, are studied with the use 
of simplest van der Waals equation of state 
(VdW-EOS). Possibility of demonstrable and 
semi- analytical description of thermo- and hy- 
drodynamics of the process is main advantage 
of this EOS. The simplest self-similar ID hy- 
drodynamics of isentropic expansion of semi- 
infinite VdW-WDM in planar geometry is an- 



alyzed. 

Caloric EOS in one-phase region. It is 

known that so-called "reduced" caloric EOS in 
form e(P, p) is the only thermodynamic EOS, 
which is necessary for description of all adia- 
batic flows: i.e. shock compression; isentropic 
expansion and compression; throttling process 
(Joule-Thompson); adiabatic expansion in vac- 
uum, istant spinodal decomposition etc. 

Thermal EOS of van der Waals (VdW) fluid 
in one-phase region is well-known (see (1)). 
In the simplest case of gas of structureless 
molecules isochoric heat capacity is constant: 
cy = const. Internal energy and entropy (with- 
out insignificant constants) in one-phase region 
can be calculated as [7]: 

(P + ap 2 )(l - pb) = pT, (1) 

e(P, p) = cyT — pa = cy —(P + ap 2 ) — pa, 

P 

(2) 

s(P,p) =cyln(P + ap 2 )+4 d ln^^, (3) 

P 

where a, b - parameters of VdW EOS and 
dp — cy + 1 = const. 

Isentrope in one-phase region follows from 
(©): 

P s (p) = const - W 2 , ( 4 ) 

where 70 = Cp/cy. 

EOS in two-phase region. For equilib- 
rium states in two-phase region EOS ([!]), Q 
and (|3| should be corrected with use of Maxwell 
rule [7 . This correction will lead to pressure be- 
coming unambiguous function of temperature 
P = P(T). From here one can derive form of 
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EOS in two-phase region. 
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Using ([5| one obtains: 



e(P,p) = l^l+g(P), 
P 



(5) 

(6) 
(7) 



where functions f(P), g(P) and <fi(P) can be 
obtained numerically. 
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where x g and x\ - value of physical quantity x 
on gaseous and liquid border of phase transition 
correspondingly. 

It should be stressed that equations (|6| and 
([7]) are universal forms of EOS in two-phase 
region for all phase transitions of 1 st order: 
they are derived without referrence to partic- 
ular EOS. 

Taking into account that {§p) v = T (§p) v 
one can derive universal dependence between 
functions g(P) and (j)(P): 
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For values of all quantities should coincide on 
border between one- and two-phase regions - 
binodal, - another expressions and mutual de- 
pendence for f(P) and g(P) can be derived 
with using ^ and EOS ([2} : 



f(P) = apip g , 
g(P)=c v T-a( Pl + p g ), 

a A P \ 



g(P)=c v T--\l 



(10) 



In two-phase region Poisson adiabat can be 
obtained numerically: 
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(11) 

Advantages of van der Waals EOS. The 

main advantage of VdW EOS is that simple ex- 
pression for f(P) and g(P) and analytical de- 
pendence (10) between them can be obtained. 



P(T) in two-phase region is sufficient to de- 
scribe state of fluid in it. 

Advantage of all EOS's with two parameters 
is that all quantities can be transformed into 
ones without dimension: reduced form of EOS 
can be used. It also can be applied to functions 
f(P) and g(P). Since there are characteristic 
values for pressure and density - P cr it and p cr n 
- one can derive characteristic values for length, 
time and, what is especially important for this 
investigation, for velocity: 
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Pcrit 



Therefore hydrodynamics of van der Waals 
fluid can be described in nondimensional form. 
Hydrodynamics in planar geometry. 

Process of isentropic expansion of semi-infinite 
layer in vacuum is considered. It is described 
by equations of hydrodynamics and EOS[8[[9]: 
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(12) 



In semi-infinite planar geometry the solution is 
self-similiar: all quantities depend on variable 
£ = x/t only [SJ [9] . Taking that into account 
one can transform system (12): 



£ = u ± c; 



u = ± 



P S (P) 



dp 



(13) 



It also means, that knowledge of dependence 



where c = y (§p ) ~ speed of sound, sign (plus 

or minus) is determined by direction of mate- 
rial flow of material. In ( 13 ) we also took into 
account the fact, that in our case flow is isen- 
tropic. So we substituted EOS for Poisson adi- 
abat P = P s (p). 

"Phase freezeout". It is clearly seen on 
Figure [l] that adiabat has a breakpoint on in- 
terception with binodal. It means, that sound 
of speed has a jump at this point. For mass 
velocity is continuous function (in other case 
the law of flow conservation would be violated), 
we also have jump of variable £ from value ^ 
- right before entering two-phase region - to 
value ^ - right after it. 

For interception point is unambiguously de- 
termined by choice of adiabat, value of sound 
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Figure 1: Interception of Poisson adiabats with 
binodal. For VdW fluid p < 3p crit 



of speed jump is determined by entropy of pro- 
cess. 

From here we can conclude, that for all values 
%k < £ < all hydro- and thermodynamical 
parameters stay constant - extensive uniform 
"boiling" layer emerges. 

Example of the problem's solution with 
clearly seen "boiling" layer is depicted on Fig- 
ure [2] Initial state and "boiling" layer for the 
example are depicted as points IS and BL cor- 
respondingly on Figure [I] 
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Figure 3: Regimes of "phase freezeout". 




Figure 2: Example of flow with "boiling" layer. 

One of the parameters of phenomenon is part 
of disturbed matter, confined in "boiling" layer 
a: 



Pb 4 



Po 



Co 



(14) 



where po and cq - density and speed of sound 
in initial state, pb - density on "boiling" layer, 
c£~ and - sound of speed before and after 
entering two-phase region. 



Figure 4: Density profiles in different regimes 
of "phase freezeout". 

For example, for the case, depicted on Figure 
U a = 73.3%. 

Parameters a and - mass velocity on 
"boiling" layer - are determined by choice of 
initial state. If entropy is already chosen, we 
can say, that they determined by value of ini- 
tial pressure. 

So we identified, that the phenomenon has 
two degrees of freedom: entropy and initial 
pressure. 

Regimes of "phase freezeout". It was 

obtained, that "boiling" layer displays different 
behaviors: direction of expansion of the layer 
changes with initial conditions. There are three 
regimes of phenomenon: 

1. c\ > > Ub - the whole layer moves in 
direction, opposite to vacuum, 

2. c£ > Ub > - the layer expands in both 
directions, 

3. Ub > c£~ > - the whole layer moves 
toward vacuum. 
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On Figure [3] initial states of boundary 
regimes — and = are depicted as 
thick lines. Dashed line is binodal, doted line 
is critical adiabat. 

On Figure [4] different regimes are illustrated. 
Initial states for all three cases are depicted as 
points 1, 2 and 3 on Figure [3] For the case 1 
a = 29.9%, for the case 2 a = 13.7%, for the 
case 3 a = 5.3%. 

Conclusion. Features of the discussed 
regimes of isentropic expansion of WDM- 
material makes it possible to discuss, at least 
formally, possibility of generation of uniform 
and extensive sample of investigated substance, 
which we are interested in, just on the boiling 
boundary of liquid-gas phase transition (includ- 
ing the region of its critical point). Size of this 
zone is determined by interception point of bin- 
odal and Poisson adiabat and initial pressure. 
It grows with growth of sound of speed jump 
at binodal, which is higher at higher densities, 
and decreases with growth of initial pressure. 
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